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ABSTRACT 

We investigate the effects of the radiatively-generated tan /3-enhanced Higgs-singlet 
Yukawa couplings on the decay T — > jAi in the NMSSM, where A\ is the lightest CP-odd 
scalar. This radiative coupling is found to dominate in the case of a highly singlet Higgs 
pseudoscalar. The branching ratio for the production of such a particle is shown to be 
within a few orders of magnitude of current experimental constraints across a significant 
region of parameter space. This represents a potentially observable signal for experiments 
at present S-factories. 
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The Next to Minimal Supersymmetric extension of the Standard Model (NMSSM) is a 
well-motivated model of electroweak symmetry breaking which resolves both the hierarchy 
problem of the Standard Model (SM) and \i problem of the Minimal Supersymmetric 
extension of the Standard Model (MSSM) in a natural way [1]. The \x parameter of the 
MSSM is replaced with an additional gauge singlet Higgs superfield S and an effective 
doublet mixing term /z e ff is generated when the singlet field acquires a vacuum expectation 
value (VEV). It has long been known that the NMSSM suffers from the formation of 
electroweak scale cosmic domain walls [2], although mechanisms to resolve this problem 
have been suggested, e.g. [3]. In the NMSSM, all parameters are naturally predicted to be 
of the order the SUSY-breaking scale Msusy- 

The Higgs sector of the NMSSM may be derived from the superpotential of the model, 
given by 

W ffiggs = \SH X H 2 + nS\ (1) 

where Hi(H 2 ) is the doublet Higgs superfield which gives masses to the down-type quarks 
and leptons (up-type quarks). The corresponding soft SUSY-breaking terms are given by 

£m gg s = AA A 5$ 1 $ 2 + kA k S 3 , (2) 

where $1^ and S are the scalar components of and S respectively. At tree level only 
two further parameters are required, the ratio of doublet VEVs tan /3 = ^ and the effective 
doublet mixing parameter /i c g = Radiative corrections due to the quarks and scalar 
quarks of the third generation must also be included in order to raise the mass of the 
SM-like Higgs H x above the LEP bound of 114 GeV. \ 

The superpotential of the NMSSM exhibits a global U(1)r symmetry which is spon- 
taneously broken when S, the scalar component of S, acquires a VEV. In addition, it is 
explicitly broken by the soft trilinear couplings A\, A K [4]. The CP-odd scalar component 
of S is therefore a pseudo-Goldstone boson of this symmetry, and is massless in the limit 
A\, A K —> 0. For small values of the trilinear couplings, the lightest pseudoscalar in the 
NMSSM spectrum can therefore naturally be very light and highly gauge singlet in nature. 
Typically this requires A\ ~ 200 GeV, A K ~ 5 GeV. Such a scenario can arise within the 
context of gauge- or gaugino-mediated SUSY breaking, where both couplings are zero at 
tree level, with non-zero A\ being radiatively generated at one loop and non-zero A K at 
two loops [5]. 

For a sufficiently light A\ boson, observation in the decay T(ls) — > jAi becomes a 



1 It is also possible to evade the LEP bound if Hi decays into the lightest pseudoscalars, with branching 
ratio B(Hi — > > 0.7 [5]. This requirement leads to a lower bound on the doublet component of Ax, 

0{\ > 0.04. Since the radiative corrections are subdominant in this region, we do not include these points 
in our results. 
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possibility^] Such a signal has previously been considered in [7,8], with the pseudoscalar 
coupling to 6-quarks only through tree level singlet-doublet mixing. It has recently been 
shown that the singlet Higgs bosons also receive a direct coupling to fermions at one loop [9] . 
Although loop suppressed, this coupling is enhanced by the ratio of Higgs doublet VEVs 
tan (3 and can become competitive with tree-level effects when this parameter is large. 

In this Letter we consider the effects of such a direct coupling on the decay T — > jAi. 
Experimental searches [10] for a light Higgs boson in T decays place a 90% confidence level 
upper bound on the branching ratio B(T — > jAi)^ 1 x 10~ 4 for a light particle m^j < 8 GeV 
decaying visibly within the detector. The upper bound rises to ~ 10 -3 for heavier particles 
due to the softness of the recoil photon and cuts placed on energy deposits in the detector 
tighten the constraints to ~ 10~ 5 for a stable or invisibly decaying A\ boson [11]. 

The branching ratio for T decays through the Wilczek mechanism [8, 12] is given by 
B (T^) _ G** J2 / _ ^ p _ (3) 



Here F ~ 1/2 includes QCD corrections [13] and B(T -> = (2.48 ± 0.06)%. The 

SM-normalised pseudoscalar coupling g^. bb is given by [9] 



P ( V2(A b ) y 

A.UU 1 1 I 

V «i J 



9A,bb 



(tan/S + A^C^ + A^-^ 



(4) 



cos j3 

with O a the 2x2 orthogonal pseudoscalar mixing matrix, such that 

Ai = 0£ 1 a + 0£ 1 a s , (5) 

where a is the would-be CP-odd scalar in the MSSM limit and as is the CP-odd singlet 
Higgs boson. In addition, A^ 2 S are the one-loop non-holomorphic Yukawa couplings of the 
states ci2.s to b quarks. At zero external momentum, they may be calculated by 

\ oa 2 ,s I 

where Ab[$i, $2, S] is a Coleman- Wienberg type functional [14] of the background Higgs 
fields which encodes radiative corrections to the b quark self-energy. Here (...) denotes 



2 In principle, this decay is also possible within other Higgs singlet extensions of the MSSM, such as the 
Minimal Non-minimal Supersymmetric extension of the Standard Model (MNSSM) [6], and Eq. ([3|) is also 
valid in this case. However, the Higgs bosons of the MNSSM obey a tree level mass sum rule, so that any 
light pseudoscalar boson is accompanied by a quasi-degenerate scalar boson. Singlet-doublet mixing in the 
scalar sector typically excludes such a scenario except in the MSSM limit of the theory A — > with /j, e s 
fixed. The radiative coupling of the singlet pseudoscalar to fermions, whose effects we consider here, will 
also vanish in this limit. 
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taking the VEV of the enclosed expression. The dominant contributions to A^ 2 S are 
due to gluino-sbottom quark and chargino-stop quark loops. In the single-Higgs-insertion 
approximation, neglecting subdominant terms proportional to the weak gauge coupling a w , 
they may be given by 



A? = -^Ma/i/^^^^-^M^M^M^/i 2 ), (7) 
Ar = -^M a ^I(^M^)-^^I(^M!y) t (8) 

where Mqja are the soft squark masses, A t is the top-squark soft trilinear coupling and 
M3 is the gluino mass. The one-loop function I(x, y, z) is given by 

xy\n(x/y) +yz\n(y/z) + xz\n(z/x) 

Vi z) = 7 w w \ • 9 

[x -y)(y -z)(x- z) 

In Fig. [1] we present results from a scan over the parameters 

< A < 0.5, < A\ < 300GeV, 
-0.5 < k < 0.5, < A K < 20GeV, 

whilst fixing tan (3 = 50 and fi c s = 120 GeV. We require a light Higgs pseudoscalar < 
9 GeV along with a lightest Higgs scalar > 114 GeV, in agreement with constraints 
from LEP II. The soft-SUSY breaking parameters which enter the calculation of A^ 2 S are 
taken to be equal at Msusy = 600 GeV. The branching ratio B(T — > jAi) is plotted against 
the non-singlet fraction of A±, described by the mixing matrix element Of x . 

The threshold corrections are independent of the tree-level coupling proportional to 
the pseudoscalar mixing, and enter the expression for bb with opposing sign. For a 
relatively large non-singlet component above few %, the threshold corrections represent a 
small suppression to the branching ratio of up to ~ 10%. In the case of a highly singlet Ai 
boson, the threshold corrections become the dominant effect, producing a branching ratio 
of the order ~ 1 x 10 -6 across a significant region of parameter space. This prediction 
is found to be generic for electroweak-scale soft SUSY-breaking terms around a TeV. At 
the intersection of these regimes, the contributions cancel giving a highly suppressed decay 
rate. 

Fig. [2] shows results from a scan over the parameter range of Eq. ( fTUI) for tan (3 = 10, 
keeping \i = 120 GeV and the common soft-SUSY breaking scale Msusy = 600 GeV. 
Both the doublet-singlet mixing and threshold correction contributions to g^ lbb are tan/3 
enhanced, such that the branching ratio at low tan (3 is smaller by 1 ~ 2 orders of magnitude 
across the full parameter space. Due to their common enhancement, the relative importance 
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Figure 1: The branching ratio B(T — > jAi) vs. the non- singlet fraction 0±± at tan/3 = 50. 
The points in green (light grey) include the one-loop threshold effects points in red 
(dark grey) neglect these corrections. Here fi c g = 120 GeV and A, k, A\,A k are scanned 
over the range given in Eq. UP) . All other soft-SUSY breaking parameters are taken to 
equal Msusy = 600 GeV. Experimental bounds are shown in dark blue (black) for a stable 
or invisibly decaying pseudoscalar and in light blue (grey) for a visibly decaying particle, 
assuming here ~ 5 GeV. The full limits are strongly dependent on the value ofrriA l and 
are less restrictive by one to two orders of magnitude for a heavy A\ boson (mA 1 > 8 GeV). 



of the two terms in Eq. (j3J) varies only slowly with tan so that for all values of tan (3 ^ 5, 
minimal branching ratios are observed for singlet-doublet mixing around few x 0.1%. 
The magnitude of the branching ratio is not found to vary strongly with Msusy or A*, 
although the available parameter space consistent with out requirements mm > 114 GeV, 
< 9 GeV decreases as \x increases, such that small values of the singlet- doublet mixing 
0^ x do not appear. 

The inclusion of threshold corrections can clearly alter the phenomenology of highly 
singlet light pseudoscalars in a dramatic way, allowing for the possibility of detectable 
T — > Ai'j decays in a new corner of parameter space. In the limit of vanishing singlet- 
doublet mixing the tree level coupling of the A\ boson to r leptons also vanishes, however an 
analogous threshold correction also contributes to the Ait + t~ coupling g^ 1T + T -, through 
a wino-stau loop. The pseudoscalar is therefore predicted to decay into t + t~ pairs with 
branching ratio of order one, for 2m r < < rax, independently of the singlet-doublet 
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Figure 2: The branching ratio B(Y — > jAi) vs. the non- singlet fraction 0^ at tan/3 = 10. 
Here fi c s = 120 Ge V and Msusy = 600 Ge V, with A, k, A\, A k scanned over the range given 
in Eq. n~U\) . Points in green (light grey) include the one-loop threshold effects A^ s ; points 
in red (dark grey) neglect these corrections. Experimental bounds are shown as for Fig. Ql 



mixing. An order-of-magnitude estimate suggests that current S-factories should be sen- 
sitive to branching ratios of the order B(T — > jAi) ^ 10 -6 , for observing such a final 
state. 

At masses above ~ 9 GeV, the A\ boson can mix with the rjb meson. This can lead 
to significant enhancement or suppression of £>(Y — > 7A1) [15]. In addition, there is a 
broadening of the A\ width, and the resonance in the energy spectrum of the recoil photon 
is less sharply peaked. There has been a suggestion to search for such a light Higgs boson 
through precision tests of lepton universality in the decays of the T [16]. Such searches 
would also be sensitive to decays in the zero-mixing limit. If the A\ boson is below the 
t + t~ threshold, the dominant decay channels are ss, gg (and hence light mesons) or photon 
pairs, since the coupling to cc is tan f3 suppressed^] This remains a favourable situation for 
the clean environment of an e + e~ collider, where these final states can be reliably measured. 

Unfortunately, despite the tremendous production rates for 6-mesons at the LHC, 
a discovery of the A± boson through this mechanism appears difficult. The final state 
consists of low-energy r jets and a photon, neither of which presents a clean signal above 

3 The tan/3 suppression would also exclude the possibility of an observable signal J/ip — > jAi for 
mA 1 < 2m c in the limit of vanishing singlet-doublet mixing. 
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background activity. An alternative production mechanism has been suggested in [17], 
which considers instead the process pp — > XiXiAi i n the limit of vanishing doublet-singlet 
mixing. The possibility for observing such a signal is strongly dependent on the masses 
and decay channels of both the lightest chargino and the Ai boson. 

If both terms contributing to b i are °f similar magnitude, typically for around 
~ 0.5% mixing, detection of the Ai boson may be extremely challenging. In this case, the 
branching fraction of T — ► A17 becomes extremely suppressed. An alternative experimental 
strategy is to look for A\ pair production from Higgs boson decays [18]. In the small singlet- 
doublet mixing scenario at large tan/3, the lightest CP-even Higgs boson Hi is highly SM- 
like, and the branching fraction B(Hi — > A^x) is conservatively bounded from above at 
around ~ 10~ 3 . Associated production of the A\ boson with a chargino pair would remain 
a possibility. 

In conclusion, we have shown that the branching ratio for production of a light Higgs 
pseudoscalar in T(ls) decays does not vanish in the absence of doublet-singlet mixing. 
We found that the decay T — » 7A1 is predicted to be observable at existing experimental 
facilities if supersymmetry is broken at the TeV scale with large tan f3. In the event of a 
cancellation between the threshold corrections and tree-level mixing contributions to the 
A\bb coupling the branching ratio may be highly suppressed even though the doublet-singlet 
mixing is still significant, and further phenomenological considerations would be needed. 
We hope to return to this issue in a future communication. 
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